The examples we will cover today involve open source text datasets that have been mostly cleaned for us. We will mostly just have to work on reformatting these data sets to get them “Tidy”.
When working with our own data this may not always be the case though. Often text data has inconsistent punctuation, spelling error, acronyms, and multiple ways of referring to the same “thing”. Several key cleaning resources can be helpful for addressing these issues:
#base functions work too!
?gsub #detects a pattern and replaces match with characters specified
?grepl #detects presence or absence of a pattern and retuns a logical vector
#tidy versions
# install.packages("stringr")
library(stringr)
?str_replace
?str_detect
A table with one-token-per-row, where a token is just a meaningful “unit” of text (e.g., word, sentence, or paragraph). This format often makes the data easier to manipulate and visualize using standard “tidy” tools.
Different from how text is typically stored in many analytic text mining approaches, which often work on a “corpus” and/or a “document-term matrix”. However we will see that transfer among the types is easy if our data is first “tidy”.
Let’s try it ourselves.
#install these packages if you don't already have them
# install.packages("janeaustenr")
# install.packages("tidytext")
# install.packages("plyr")
# install.packages("dplyr")
# install.packages("magrittr")
library(tidytext)
library(janeaustenr)
library(plyr)
library(dplyr)
library(magrittr)
#browse the related vignette or others
vignette("tidytext",package="tidytext")
?austen_books
original_books <- austen_books() %>%
group_by(book) %>% #grouping df by "book"
mutate(linenumber = row_number(), #adding row numbers as a new variable
chapter = cumsum(str_detect(text, regex("^chapter [\\divxlc]",
ignore_case = TRUE)))) %>%
ungroup()
# the above is creating a new chapter variable by searching the "text" character variable for the string "chapter" and Roman Numeral characters (as specified in regex language here "[\\divxlc]")
head(original_books)
tail(original_books)
slice(original_books, 1:500) %>% View
#restructuring as a one-token-per-row format. unnest() does that for us
tidy_books <- original_books %>%
unnest_tokens(word, text)
#The default tokenizing is for words, but other options include characters, ngrams, sentences, lines, paragraphs, or separation around a regex pattern.
data("stop_words")
head(stop_words)
cleaned_books <- tidy_books %>%
anti_join(stop_words)
#most common words in all the books?
cleaned_books %>%
count(word, sort = TRUE)
#Examining most common words across all books and by book:
cleaned_books_2 <- cleaned_books %>%
mutate(word = str_extract(word, "[[:alpha:]]+")) %>%
count(word,sort = TRUE) %>% #counting the number of occurrences by group (word in this case)
left_join(cleaned_books,.) %>%
rename(n_allbooks = n) %>%
group_by(book,word) %>%
mutate(n_bybook = n()) %>%
arrange(book, desc(n_bybook)) %>%
ungroup()
##str_extract = extract strings that match the patten spcified, in this case, extracting all alphabetic characters regardless of capitalization. So Numbers, symbols (i.e.,"_" and "*" often used to flag italics or bold will be removed-- we don't want don't want "Angry" to appear as a diff word than it's italicized form "_Angry_")
#Chart of top 5 words across all books
#install.packages("ggplot2")
library(ggplot2)
#doing all in one tidy flow
cleaned_books %>%
count(word, sort = TRUE) %>%
filter(n > 600) %>%
mutate(word = reorder(word, n)) %>%
ggplot(aes(word, n)) +
geom_col() +
xlab(NULL) +
coord_flip()
#using new var we created above
#no scientifict notation
options(scipen = 10000)
cleaned_books_2 %>%
filter(n_allbooks > 600, !duplicated(word)) %>%
mutate(word = reorder(word, n_allbooks)) %>%
ggplot(aes(word, n_allbooks)) +
geom_col() +
xlab(NULL) +
coord_flip()
#What about the top 5 words only in each book?
cleaned_books_2 %>%
group_by(book) %>%
filter(!duplicated(word)) %>%
arrange(desc(n_bybook)) %>%
top_n(3, n_bybook) %>%
ungroup %>%
arrange(book, n_bybook) %>%
ggplot(aes(word, n_bybook)) +
geom_col() +
xlab(NULL) +
coord_flip() +
facet_wrap(~book) +
ylab('# of Times Used') +
xlab('word')
#other ways to visualize??
We can use the tools of text mining to assess the emotional content of text:
Flowchart: using tidytext for sentiment analysis
#Four sentiment lexicons are in the tidytext package in the sentiments dataset.
head(sentiments)
## # A tibble: 6 x 4
## word sentiment lexicon score
## <chr> <chr> <chr> <int>
## 1 abacus trust nrc NA
## 2 abandon fear nrc NA
## 3 abandon negative nrc NA
## 4 abandon sadness nrc NA
## 5 abandoned anger nrc NA
## 6 abandoned fear nrc NA
table(sentiments$lexicon)
##
## AFINN bing loughran nrc
## 2476 6788 4149 13901
#How many distinct "emotions" are words categorized into?
table(sentiments$sentiment)
##
## anger anticipation constraining disgust fear
## 1247 839 184 1058 1476
## joy litigious negative positive sadness
## 689 903 10461 4672 1191
## superfluous surprise trust uncertainty
## 56 534 1231 297
For more on the sentiment lexicons contained in this dataset and how they differ checkout some details here
Keep in mind that these methods do not take into account qualifiers before a word, such as in “no good” or “not true”; a lexicon-based method like this is based on unigrams only. For many kinds of text (like the narrative examples below), there are not sustained sections of sarcasm or negated text, so this is not an important effect.
Also, the size of the chunk of text that we use to add up sentiment scores can effect analysis. A text the size of many paragraphs can often have positive and negative sentiment averaged out to about zero, while sentence-sized or paragraph-sized text often works better.
Let’s try some sentiment analysis ourselves
#Let’s look at the words with a joy score from the NRC lexicon. What are the most common joy words in Emma?
nrcjoy <- get_sentiments("nrc") %>%
filter(sentiment == "joy")
head(nrcjoy)
## # A tibble: 6 x 2
## word sentiment
## <chr> <chr>
## 1 absolution joy
## 2 abundance joy
## 3 abundant joy
## 4 accolade joy
## 5 accompaniment joy
## 6 accomplish joy
#"get_sentiments" is just a helper fxn that retrives a specific lexicon with only relevant columns
# equivalent too: filter(sentiments, lexicon == 'nrc' & sentiment == 'joy') %>% select(word, sentiment)
emma_joy <- tidy_books %>%
filter(book == "Emma") %>%
semi_join(nrcjoy) %>%
count(word, sort = TRUE)
head(emma_joy)
## # A tibble: 6 x 2
## word n
## <chr> <int>
## 1 good 359
## 2 young 192
## 3 friend 166
## 4 hope 143
## 5 happy 125
## 6 love 117
#notice how we used the original tidy books dataset as opposed to to "cleaned_books_2" which we later made. Why? Stop words (i.e., "good", "young") are often overlapping with sentient phrases.
#We could also examine how sentiment changes during each novel.
#Let’s find a sentiment score for each word using the Bing lexicon, then count the number of positive and negative words in defined sections of each novel.
# a little on bing :
#categorizes words in a binary fashion to be "postiive" or "negative"
bing <- get_sentiments("bing")
head(bing)
## # A tibble: 6 x 2
## word sentiment
## <chr> <chr>
## 1 2-faced negative
## 2 2-faces negative
## 3 a+ positive
## 4 abnormal negative
## 5 abolish negative
## 6 abominable negative
#install.packages("tidyr")
library(tidyr) #we will use this for reformatting df (see ?spread)
janeaustensentiment <- tidy_books %>%
inner_join(get_sentiments("bing")) %>% #return tidy_book rows with matching values in bing sentiment df
count(book, index = linenumber %/% 100, sentiment) %>% #counting number of positive/neg lines wihtin each index chunk
spread(sentiment, n, fill = 0) %>%
mutate(sentiment = positive - negative) #simple measure for each index as # positive words - # of neg words.
#indexing will be used to parse up the novel and examine sentiment change throughout
#index 0 represents the first 100 lines, index 1 represents the 2nd 100 lines, etc.
#"%/%" does "integer" division for us
5/2
## [1] 2.5
#vs
5%/%2
## [1] 2
#let's look at sentiment in the different novels
ggplot(data = janeaustensentiment, aes(x = index, y = sentiment, fill = book)) +
geom_bar(alpha = 0.8, stat = "identity", show.legend = FALSE) +
facet_wrap(facets = ~ book, ncol = 2, scales = "free_x")
#We could also examine the most common positive & negative words in the Jane Austen novels
bing_word_counts <- tidy_books %>%
inner_join(get_sentiments("bing")) %>%
count(word, sentiment, sort = T) %>%
ungroup()
#format for graphing
tmp <- bing_word_counts %>%
filter(n > 150) %>%
mutate(n = ifelse(sentiment == "negative", -n, n)) %>% #just making counts for "neg" words negative for graphing
mutate(word = reorder(word, n))
ggplot(data = tmp, mapping = aes(x = word, y = n, fill = sentiment)) +
geom_bar(alpha = 0.8, stat = "identity") +
labs(y = "Contribution to sentiment", x = NULL) +
coord_flip()
#This lets us spot an anomaly in the sentiment analysis; the word “miss” is coded as negative but it is used as a title for young, unmarried women in Jane Austen’s works. If it were appropriate for our purposes, we could easily add “miss” to a custom stop-words list using bind_rows(). We could implement that with a strategy such as this.
custom_stop_words <- bind_rows(data_frame(word = c("miss"),
lexicon = c("custom")),
stop_words)
custom_stop_words
## # A tibble: 1,150 x 2
## word lexicon
## <chr> <chr>
## 1 miss custom
## 2 a SMART
## 3 a's SMART
## 4 able SMART
## 5 about SMART
## 6 above SMART
## 7 according SMART
## 8 accordingly SMART
## 9 across SMART
## 10 actually SMART
## # ... with 1,140 more rows
So far we’ve analyzed text in the tidy format (one-token-per-doc-per-row). Many useful descriptive analyses can be done from this format, but many NLP analyses require a different format:
Flowchart: Tidy Text to other Formats
Let’s examine one of the most common non-tidy structures text mining packages use: the document-term matrix (DTM). This is a matrix in which:
Let’s convert our tidy_books df to a dtm so we can impliment one common NLP model: Latent Dirichlet allocation (i.e., LDA or topic modeling)
#install the tm and topicmodels package if you don't have them
#install.packages("tm")
#install.packages("topicmodels")
library(tm)
library(topicmodels)
austen_dtm <- austen_books() %>%
unnest_tokens(word, text) %>%
count(book, word) %>%
cast_dtm(book, word, n)
austen_dtm <- austen_books() %>%
group_by(book) %>%
ungroup() %>%
unnest_tokens(word, text) %>%
anti_join(stop_words) %>%
mutate(word = str_extract(word, "[[:alpha:]]+")) %>% #taking out chp numbers (and others!)
count(book, word) %>%
cast_dtm(book, word, n)
#let's take a look at the cast_dtm fxn
?cast_dtm
austen_dtm
## <<DocumentTermMatrix (documents: 6, terms: 13311)>>
## Non-/sparse entries: 36350/43516
## Sparsity : 54%
## Maximal term length: NA
## Weighting : term frequency (tf)
Topic modeling helps us devide collections of documents into groups so we can understand their content better.
Topic modeling is:
LDA is a popular topic modeling method that treats:
Suppose a vandal has broken into your study and torn apart four of your books:
This vandal has torn the books into individual chapters, and left them in one large pile. How can we restore these disorganized chapters to their original books? This is a challenging problem since the individual chapters are unlabeled: we don’t know what words might distinguish them into groups. We’ll thus use topic modeling to discover how chapters cluster into distinct topics, each of them (presumably) representing one of the books.
We’ll retrieve the text of these four books using the gutenbergr package
titles <- c("Twenty Thousand Leagues under the Sea", "The War of the Worlds",
"Pride and Prejudice", "Great Expectations")
#install.packages("gutenbergr")
library(gutenbergr)
books <- gutenberg_works(title %in% titles) %>%
gutenberg_download(meta_fields = "title")
# divide into documents, each representing one chapter
by_chapter <- books %>%
group_by(title) %>%
mutate(chapter = cumsum(str_detect(text, regex("^chapter ", ignore_case = TRUE)))) %>%
ungroup() %>%
filter(chapter > 0) %>%
unite(document, title, chapter)
# split into words
by_chapter_word <- by_chapter %>%
unnest_tokens(word, text)
# find document-word counts
word_counts <- by_chapter_word %>%
anti_join(stop_words) %>%
count(document, word, sort = TRUE) %>%
ungroup()
word_counts
## # A tibble: 104,721 x 3
## document word n
## <chr> <chr> <int>
## 1 Great Expectations_57 joe 88
## 2 Great Expectations_7 joe 70
## 3 Great Expectations_17 biddy 63
## 4 Great Expectations_27 joe 58
## 5 Great Expectations_38 estella 58
## 6 Great Expectations_2 joe 56
## 7 Great Expectations_23 pocket 53
## 8 Great Expectations_15 joe 50
## 9 Great Expectations_18 joe 50
## 10 The War of the Worlds_16 brother 50
## # ... with 104,711 more rows
chapters_dtm <- word_counts %>%
cast_dtm(document, word, n)
chapters_dtm
## <<DocumentTermMatrix (documents: 193, terms: 18215)>>
## Non-/sparse entries: 104721/3410774
## Sparsity : 97%
## Maximal term length: 19
## Weighting : term frequency (tf)
chapters_lda <- LDA(chapters_dtm, k = 4, control = list(seed = 1234))
chapters_lda
## A LDA_VEM topic model with 4 topics.
#setting the seed is important to make sure our results are reproducible
#Let's extract the per-topic-per-word probabilities.
chapter_topics <- tidy(chapters_lda, matrix = "beta")
chapter_topics
## # A tibble: 72,860 x 3
## topic term
## <int> <chr>
## 1 1 joe
## 2 2 joe
## 3 3 joe
## 4 4 joe
## 5 1 biddy
## 6 2 biddy
## 7 3 biddy
## 8 4 biddy
## 9 1 estella
## 10 2 estella
## # ... with 72,850 more rows, and 1 more variables: beta <dbl>
#turned the model into a one-topic-per-term-per-row format. For each combination, the model computes the probability of that term being generated from that topic. For example, the term “joe” has an almost zero probability of being generated from topics 1, 2, or 3, but it makes up 1.45% of topic 4
#what are the top 5 terms within each topic?
top_terms <- chapter_topics %>%
group_by(topic) %>%
top_n(5, beta) %>%
ungroup() %>%
arrange(topic, -beta)
top_terms
## # A tibble: 20 x 3
## topic term beta
## <int> <chr> <dbl>
## 1 1 elizabeth 0.014107538
## 2 1 darcy 0.008814258
## 3 1 miss 0.008706741
## 4 1 bennet 0.006947431
## 5 1 jane 0.006497512
## 6 2 captain 0.015507696
## 7 2 nautilus 0.013050048
## 8 2 sea 0.008850073
## 9 2 nemo 0.008708397
## 10 2 ned 0.008030799
## 11 3 people 0.006797400
## 12 3 martians 0.006512569
## 13 3 time 0.005347115
## 14 3 black 0.005278302
## 15 3 night 0.004483143
## 16 4 joe 0.014450300
## 17 4 time 0.006847574
## 18 4 pip 0.006817363
## 19 4 looked 0.006365257
## 20 4 miss 0.006228387
#let's visualize
library(ggplot2)
top_terms %>%
mutate(term = reorder(term, beta)) %>%
ggplot(aes(term, beta, fill = factor(topic))) +
geom_col(show.legend = FALSE) +
facet_wrap(~ topic, scales = "free") +
coord_flip()
These topics are pretty clearly associated with the four books! There’s no question that the topic of “captain”, “nautilus”, “sea”, and “nemo” belongs to Twenty Thousand Leagues Under the Sea, and that “jane”, “darcy”, and “elizabeth” belongs to Pride and Prejudice. We see “pip” and “joe” from Great Expectations and “martians”, “black”, and “night” from The War of the Worlds. We also notice that, in line with LDA being a “fuzzy clustering” method, there can be words in common between multiple topics, such as “miss” in topics 1 and 4, and “time” in topics 3 and 4.
#Let's extract the document-topic-probabilities
#These may help us to see which topics are associated with each document
chapters_gamma <- tidy(chapters_lda, matrix = "gamma")
chapters_gamma
## # A tibble: 772 x 3
## document topic gamma
## <chr> <int> <dbl>
## 1 Great Expectations_57 1 0.00001351886
## 2 Great Expectations_7 1 0.00001470726
## 3 Great Expectations_17 1 0.00002117127
## 4 Great Expectations_27 1 0.00001919746
## 5 Great Expectations_38 1 0.35444031280
## 6 Great Expectations_2 1 0.00001723723
## 7 Great Expectations_23 1 0.55072412419
## 8 Great Expectations_15 1 0.01682502911
## 9 Great Expectations_18 1 0.00001272044
## 10 The War of the Worlds_16 1 0.00001084337
## # ... with 762 more rows
#Each of these values is an estimated proportion of words from that document that are generated from that topic. For example, the model estimates that each word in the Great Expectations_57 document has only a 0.00135% probability of coming from topic 1 (Pride and Prejudice).
#Now that we have these topic probabilities, we can see how well our unsupervised learning did at distinguishing the four books. We’d expect that chapters within a book would be found to be mostly (or entirely), generated from the corresponding topic.
#First we re-separate the document name into title and chapter, after which we can visualize the per-document-per-topic probability for each
chapters_gamma <- chapters_gamma %>%
separate(document, c("title", "chapter"), sep = "_", convert = TRUE)
chapters_gamma %>%
mutate(title = reorder(title, gamma * topic)) %>%
ggplot(aes(factor(topic), gamma)) +
geom_boxplot() +
facet_wrap(~ title)
#We notice that almost all of the chapters from Pride and Prejudice, War of the Worlds, and Twenty Thousand Leagues Under the Sea were uniquely identified as a single topic each.
#It does look like some chapters from Great Expectations were somewhat associated with other topics.
#Are there any cases where the topic most associated with a chapter belonged to another book? First we’d find the topic that was most associated with each chapter using top_n(), which is effectively the “classification” of that chapter.
chapter_classifications <- chapters_gamma %>%
group_by(title, chapter) %>%
top_n(1, gamma) %>%
ungroup()
chapter_classifications
## # A tibble: 193 x 4
## title chapter topic gamma
## <chr> <int> <int> <dbl>
## 1 Great Expectations 23 1 0.5507241
## 2 Pride and Prejudice 43 1 0.9999610
## 3 Pride and Prejudice 18 1 0.9999654
## 4 Pride and Prejudice 45 1 0.9999038
## 5 Pride and Prejudice 16 1 0.9999466
## 6 Pride and Prejudice 29 1 0.9999300
## 7 Pride and Prejudice 10 1 0.9999203
## 8 Pride and Prejudice 8 1 0.9999134
## 9 Pride and Prejudice 56 1 0.9999337
## 10 Pride and Prejudice 47 1 0.9999506
## # ... with 183 more rows
#We can then compare each to the “consensus” topic for each book (the most common topic among its chapters), and see which were most often misidentified.
book_topics <- chapter_classifications %>%
count(title, topic) %>%
group_by(title) %>%
top_n(1, n) %>%
ungroup() %>%
transmute(consensus = title, topic)
chapter_classifications %>%
inner_join(book_topics, by = "topic") %>%
filter(title != consensus)
## # A tibble: 2 x 5
## title chapter topic gamma consensus
## <chr> <int> <int> <dbl> <chr>
## 1 Great Expectations 23 1 0.5507241 Pride and Prejudice
## 2 Great Expectations 54 3 0.4803234 The War of the Worlds
#We see that only two chapters from Great Expectations were misclassified, as LDA described one as coming from the “Pride and Prejudice” topic (topic 1) and one from The War of the Worlds (topic 3).
Some additional resources to dive into later on your text mining journey:
#Document topic probabilities can be used as predictors
chapters_gamma <- tidy(chapters_lda, matrix = "gamma") %>%
separate(document, c("title", "chapter"), sep = "_", convert = TRUE) %>%
mutate(topic = paste0("topic","_",topic))
chps_topic_probability <- spread(chapters_gamma, topic, gamma) %>%
mutate(book = as.factor(str_replace_all(tolower(title)," ","_")))
#install.packages(MASS)
library(MASS)
# find the linear combinations of the original variables (in this case doc-topic-probabilities) that gives the best possible separation between the groups (i.e., books as the data is currently formatted)
chps_topic_probability <- chps_topic_probability %>%
mutate(id = row.names(.))
train <- chps_topic_probability %>%
sample_frac(.80)
test <- chps_topic_probability %>%
filter(!id %in% train$id)
lda.fit <- lda(book ~ topic_1 + topic_2 + topic_3 + topic_4 ,
data = train)
#note lda fxn here performs linear discriminant analysis-- not LDA Latent Dirichlet Allocation
lda.pred=predict(lda.fit, test)
names(lda.pred)
## [1] "class" "posterior" "x"
lda.class=lda.pred$class
table(lda.class,test$book)
##
## lda.class great_expectations
## great_expectations 15
## pride_and_prejudice 1
## the_war_of_the_worlds 0
## twenty_thousand_leagues_under_the_sea 0
##
## lda.class pride_and_prejudice
## great_expectations 0
## pride_and_prejudice 11
## the_war_of_the_worlds 0
## twenty_thousand_leagues_under_the_sea 0
##
## lda.class the_war_of_the_worlds
## great_expectations 0
## pride_and_prejudice 0
## the_war_of_the_worlds 5
## twenty_thousand_leagues_under_the_sea 0
##
## lda.class twenty_thousand_leagues_under_the_sea
## great_expectations 0
## pride_and_prejudice 0
## the_war_of_the_worlds 0
## twenty_thousand_leagues_under_the_sea 7
lda.class[1:20]
## [1] great_expectations great_expectations great_expectations
## [4] great_expectations great_expectations great_expectations
## [7] pride_and_prejudice great_expectations great_expectations
## [10] great_expectations great_expectations great_expectations
## [13] great_expectations great_expectations great_expectations
## [16] great_expectations pride_and_prejudice pride_and_prejudice
## [19] pride_and_prejudice pride_and_prejudice
## 4 Levels: great_expectations ... twenty_thousand_leagues_under_the_sea
test[1:20,'book']
## # A tibble: 20 x 1
## book
## <fctr>
## 1 great_expectations
## 2 great_expectations
## 3 great_expectations
## 4 great_expectations
## 5 great_expectations
## 6 great_expectations
## 7 great_expectations
## 8 great_expectations
## 9 great_expectations
## 10 great_expectations
## 11 great_expectations
## 12 great_expectations
## 13 great_expectations
## 14 great_expectations
## 15 great_expectations
## 16 great_expectations
## 17 pride_and_prejudice
## 18 pride_and_prejudice
## 19 pride_and_prejudice
## 20 pride_and_prejudice
mean(lda.class==test$book)
## [1] 0.974359
##Handy Function for Topic Visualization
#' @title Interactive Tool for visualizing Topics and Per-Word-pertopic-probabilities
#' @param file Fitted LDA output, corpus file: row for each token-document combo with count of token within document ,document term matrix
#' @return Shiny visual of all topics and associated words
topicmodels_json_ldavis <- function(fitted, corpus, doc_term){
# Required packages
library(topicmodels)
library(dplyr)
library(stringi)
library(tm)
library(LDAvis)
# Find required quantities
phi <- posterior(fitted)$terms %>% as.matrix
theta <- posterior(fitted)$topics %>% as.matrix
vocab <- colnames(phi)
doc_length <- vector()
for (i in 1:length(corpus)) {
temp <- paste(corpus[[i]]$content, collapse = ' ')
doc_length <- c(doc_length, stri_count(temp, regex = '\\S+'))
}
#temp_frequency <- inspect(doc_term)
freq_matrix <- data.frame(ST = colnames(doc_term),
Freq = colSums(as.matrix(doc_term)))
rm(temp_frequency)
# Convert to json
json_lda <- LDAvis::createJSON(phi = phi, theta = theta,
vocab = vocab,
doc.length = doc_length,
term.frequency = freq_matrix$Freq)
return(json_lda)
}
topics_json_lda <- topicmodels_json_ldavis(chapters_lda,
tm::Corpus(tm::VectorSource(distinct(word_counts, document)$document)),
chapters_dtm)
serVis(topics_json_lda, open.browser = TRUE)
We can take advantage of parrallel processing in R to find the optimal number of topics
#Selecting the optimal number of topics using 5-fold cross-validation with parallel processing
#Caution- this will take a few mins
#install the below packages before running if you do not already have them
library(parallel)
library(foreach)
library(doParallel)
cluster <- makeCluster(detectCores(logical = TRUE) - 1) # leave at least one CPU spare...
registerDoParallel(cluster)
clusterEvalQ(cluster, {
library(topicmodels)
})
# set parameters for LDA below
folds <- 5
candidate_k <- c(2, 3, 4, 5, 6, 7, 10)
burnin = 1000
iter = 1000
keep = 50
n <- nrow(chapters_dtm) # doc term matrix created above
splitfolds <- sample(1:folds, n, replace = TRUE)
clusterExport(cluster, c("chapters_dtm", "burnin", "iter", "keep", "splitfolds", "folds", "candidate_k", "n"))
system.time({
results <- foreach(j = 1:length(candidate_k), .combine = rbind) %dopar%{
k <- candidate_k[j]
results_1k <- matrix(0, nrow = folds, ncol = 2)
colnames(results_1k) <- c("k", "perplexity")
for(i in 1:folds){
train_set <- chapters_dtm[splitfolds != i , ] #replace with approprite dtm
valid_set <- chapters_dtm[splitfolds == i, ]
fitted <- LDA(train_set, k = k, method = "Gibbs",
control = list(burnin = burnin, iter = iter, keep = keep, seed = 1234) )
results_1k[i,] <- c(k, perplexity(fitted, newdata = valid_set))
}
return(results_1k)
}
})
stopCluster(cluster)
results_df <- as.data.frame(results)
#Selecting optimal number of topics given lowest perplexity
ggplot(results_df, aes(x = k, y = perplexity)) +
geom_point() +
geom_smooth(se = FALSE) +
ggtitle("5-fold cross-validation of topic modelling with the 'Chapters' data",
"(ie five different models fit for each candidate number of topics)") +
labs(x = "Candidate number of topics", y = "Perplexity when fit to the hold-out set")
#if you run you'll see 4/5 is the optimal number of topics!